Supra-oscillatory critical temperature dependence of Nb-Ho bilayers 
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We investigate the critical temperature T c of a thin s-wave superconductor (Nb) proximity coupled to a helical 
rare earth ferromagnet (Ho). As a function of the Ho layer thickness, we observe multiple oscillations of T c 
superimposed on a slow decay, that we attribute to the influence of the Ho on the Nb proximity effect. Because 
of Ho inhomogeneous magnetization, singlet and triplet pair correlations are present in the bilayers. We take 
both into consideration when solving the self consistent Bogoliubov-de Gennes equations, and we observe a 
reasonable agreement. We also observe non-trivial transitions into the superconducting state, the zero resistance 
state being attained after two successive transitions which appear to be associated with the magnetic structure 
of Ho. 



PACS numbers: 74.45+c, 74.62.-c,74.20.Rp,74.78.Fk 

Introduction - In recent years, the interest in 
superconductor-ferromagnet (S-F) hybrids has greatly 
intensified, partly due to a number of groups reporting the 
existence of long-range spin polarized supercurrents [1-5]. 
It has been known for many years that superconducting 
correlations can penetrate from a superconductor into a 
non- superconducting non-magnetic metal (N) via Andreev 
reflection at the S-N interface, where Cooper pairs penetrate 
the normal metal N as coherent electron-hole Andreev pairs. 
As long as the phase coherence between the electron and the 
hole is preserved, the normal metal will exhibit supercon- 
ducting features, such as the ability to carry a dissipationless 
supercurrent. At low temperatures (T ~ IK), such a 
supercurrent can penetrate as far as 1 /im into a clean normal 
metal. In a S-F junction, the penetration of superconducting 
correlations into the F is much shorter ranged (typically less 
than 10 nm even for weak F's). Two mechanisms account 
for this: the dephasing effect of the vector potential due to 
the internal magnetic field (orbital effect); and the exchange 
interaction, which acts to impose a parallel alignment of 
the spins in the Andreev pairs. This parallel alignment is at 
odds with the Andreev pairs' singlet state, which requires 
antiparallel spins. Due to the energy splitting between the 
spin-up and spin-down conduction bands in the F, the singlet 
pair amplitude oscillates in position within the F. This induces 
oscillations in the supercurrent of a S-F-S Josephson junction 
and in the critical temperature (T c ) of a S-F bilayer. In 
experiments, many oscillation periods have been observed 
in the supercurrent of S-F-S junctions [6-8]; however, S-F 
bilayers only exhibited reentrant superconductivity followed 
by one oscillation and then saturation [9-13]. In this Letter 
we report multiple T c oscillations with dp in Nb-Ho systems, 
which cannot be explained within the standard singlet pairing 
framework and that we attribute to the influence of the 
inhomogeneous Ho magnetisation on the Nb proximity effect. 

It was recently found that spin-polarized triplet Cooper 
pairs, in which the spins align parallel rather than antiparallel, 



can be generated by coupling an s-wave superconductor to an 
inhomogeneous F [14-16]. This striking possibility opens up 
a whole new field of research, in which superconductivity and 
ferromagnetism act cooperatively to create long-range, spin 
polarized, dissipation-free currents. Of most relevance to this 
paper, the intrinsic conical magnetism of Ho has been shown 
to generate a triplet supercurrent which was able to penetrate 
a 16-nm- thick film of Co [1]. Here we report experiments to 
examine the influence of the triplet pair correlations on the 
proximity effect of a thin s-wave superconductor. We mea- 
sure the critical temperature T c of Nb-Ho thin film bilayers 
whilst varying the Ho layer thickness dn m the 1.1 — 15 nm 
range. We observe multiple T c vs dn oscillations superim- 
posed on a decaying background. These arise from both the 
spin-splitting effect of the exchange field and from its inhomo- 
geneity, responsible for the spin-polarized triplet correlations. 
Furthermore, the zero resistance state in the R(T) curves is at- 
tained after two separate sharp transitions for certain Ho layer 
thicknesses. This suggests an extreme sensitivity of the prox- 
imity effect to the details of the Ho magnetic structure. 

Samples and methods - Holmium is a rare-earth ferro- 
magnet, which possesses a helical magnetic structure below 
133 K. At low temperature (T < 19 K), Ho is a spiral 
conical ferromagnet: its atomic moment of 10.34 /i^ forms 
an angle a = 80° with the crystallographic c-axis which 
is perpendicular to the sample plane, thus having a small 
constant out-of-plane component. Its in-plane component 
rotates along the c-axis, making a helix with an average turn 
angle of fj = 30° per atomic layer, giving a helical period of 
A = 3.34 nm (see fig. 2, inset). Thin poly crystalline films of 
cold- sputtered Ho present a primary (002) phase, the c-axis 
being normal to the sample surface, and a secondary (101) 
phase [18]. The Ho helical behaviour is thus relevant for the 
proximity effect even in poly crystalline films, and is clearly 
required for the generation of triplet supercurrents in SFS 
junctions [1]. The magnetisation loops show a behaviour 
qualitatively similar to that of a bulk Ho sample, but with a 
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FIG. 1. (Color online) R(T) curves with different in-plane fields 
applied B = 0,1,2,3,4,5T (from right to left, see arrow) for 
(Iho — 4.4 nm (a), cIho — 5.4 nm (b), (Iho — 7.4 nm (c) of Tl, 
and for a control Nb sample (d). 



lower saturation moment, suggesting that significant mag- 
netic inhomogeneity is present in sputtered poly crystalline 
films [17]. Two series of Nb-Ho-Cu trilayer thin films were 
prepared (labeled here as Tl and T2). The Nb and Cu 
thicknesses were kept constant, while the Ho layer thickness 
was varied between samples. The Tl and T2 structures are: 
Nb(16)-Ho[4.4-15]-Cu(20) and Nb(15)-Ho[l.l-13]-Cu(20), 
where the layer thicknesses are noted in nm units in the 
parentheses and are known with an incertitude of ±0.3 nm. 
The Cu layer was included to prevent oxidation of the Ho. 
The trilayers were deposited at ambient temperature by 
d.c. magnetron sputtering in 1.5 x 10 ~ 2 mbar of Ar onto 
oxidized Si substrates. The base pressure of the chamber 
was 10 _9 mbar, with an O2 partial pressure < 10 -11 mbar. 
The samples were mounted on a rotating sample holder, 
whose angular velocity below the magnetron determines the 
thickness of the deposited layer: all the samples in each 
series were therefore deposited without breaking the vacuum 
[19]. While the growth conditions of all the samples in 
one run are exactly the same, the only difference being the 
monotonically increasing Ho thickness, variations between 
different runs can be expected due to changes in the residual 
oxygen pressure. The film resistance (R) as a function of 
temperature (T) was measured in a four-point configuration 
with an alternating current / = ±50/iA. The R(T) curves at 
zero magnetic field were measured in a liquid 4 He dip probe 
down to T = 4.2 K, while the R(T) curves in the field range 
— 5T were measured in a pumped 4 He cryostat whose 
minimum temperature was 1.7 K. 

Results - The T c values are extracted from the R(T) 
curves; representative data are shown in figs. 1(a)- 1(c) for 



three samples with in-plane magnetic field applied in the 
0-5 T range. Two different cases are observed when varying 
the Ho layer thickness: either the transition is sharp with a 
width below 0.4 K [fig. 1(a)], or two transitions are present 
with a separation as large as 1.3 K at B = [fig. 1(c)]. To 
exclude the possibility that the double transition is due to 
inhomogeneity in the Nb layer thickness, we have deposited 
a 15-nm-thick control sample of Nb. Its superconducting 
transition [fig. 1(d)] is sharp with a transition width of 
< 0.15 K. Control Ho films, deposited over thin (4nm) Nb 
layers also show a monotonic R(T) curve in the experimental 
temperature range [17]. The complex transitions observed 
are therefore due to the interaction of the superconducting 
layer with the magnetic Ho layer: the films are structurally 
homogeneous, the heterogeneous electrical properties be- 
ing seemingly generated by the heterogeneous magnetic 
structure. To follow the dependence of the two transitions 
observed, we record the temperatures corresponding to 98%, 
50% and 10% of the normal state resistance (Rn) measured 
at 8.5 K (labeled T C5 98, T C5 5o and T C5 io respectively). T C5 io 
and T C5 5o have similar values, and correspond to the stronger 
proximity suppression of T c , leading to the zero-resistance 
superconducting state. T C5 5o is the critical temperature of the 
transition responsible for more than 80% of the resistance 
drop and in many samples the only transition present in the 
R(T) curves. Moreover, as we will show in fig. 4, T C)5 o is 
strongly affected by the magnetic field 'unzipping' the Ho 
spiral magnetisation. In fig. 2 we show T c 50 as a function of 




FIG. 2. (Color online) Superconducting critical temperature T c ,50 vs 
Ho thickness duo of the Tl series (red, d s = 16nm) and the T2 se- 
ries (blue, d s = 15 nm). The oscillations are superimposed on an ex- 
ponential decay (dotted lines). The error on T C) 5o is below ±50 mK 
(smaller than the dot), while the error on the Ho layer thickness is 
±0.3 nm; the error bar, which is the same for all the points, is shown 
as an example for the first point in each series. The transition for 
the T2 sample with dHo = 9.6 nm is not complete (T c ,5o < 4.2 K); 
however, the beginning of the transition being visible, we can esti- 
mate 3K< T c , 50 < 4.2 K. Inset: Helical magnetic structure of bulk 
Ho. 
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FIG. 3. (Color online) Solution of the BdG equations for the Tl 
series (red, d s = 16 nm) and the T2 series (blue, d s = 15 nm). 
The theoretical curve for Tl was calculated using £ho = 3.2 nm 
and an overall T c shift of 1.4 K; the curve for T2 was calculated 
using £ho = 2.1 nm and an overall T c shift of 2 K. For both curves 
T r ° = 9 K was used. 



duo for Tl and T2. In both series T C)50 oscillates, with a pe- 
riod Adno — 1-8 ± 0.5 nm for Tl and Adno = 1.5 =b 0.5 nm 
for T2. The amplitude of the oscillations can be as large 
as 0.5 K in Tl and 1 K in T2, and is much larger than the 
estimated error in the T c measurement (indicated on the 
left-most data point in each case). The oscillations are 
superimposed on an exponential decay exp(— da oft) with 
characteristic length scales I = 47 nm for Tl and I = 33 nm 
for T2. 

Oscillatory T c - In our Nb-Ho bilayers we can observe mul- 
tiple oscillations, superimposed on a slow decay, which could 
correspond to the beginning of a reentrance similar to those 
observed in homogeneous F-S films (see also the larger range 
data presented in [17]). The clearer oscillations in T2 are pos- 
sibly induced by the thinner Nb layer, on which the Ho has 
a stronger effect, leading to smaller T c and larger oscillation 
amplitudes [9]. 

To compare our results with theory, we have calculated 
the transition temperature of Nb-Ho bilayers using a fully 
microscopic procedure, based on the exact numerical di- 
agonalization of the self-consistent Bogoliubov-de Gennes 
(BdG) equations [20-22]. This method can effectively handle 
the geometrical effects inherent to finite systems, as well 
as the triplet correlations generated by the presence of the 
nonuniform magnetization in Ho [23]. The m = and 
m = ±1 triplet amplitudes are all found to be nonvanishing, 
and similar in importance. The bilayers are theoretically 
modeled as being infinite in the transverse direction. The 
helical magnetic structure of Ho is described by a spiral local 
exchange field whose angle with the c-axis a, turning angle 



j3 and helical period A are given in the literature. In general, 
whenever possible, we set the relevant parameters to their 
corresponding experimental values, keeping the adjustable 
ones to a minimum. Thus, we set the Ho and Nb layers 
thicknesses duo and d s , and the critical temperature T® of 
a bulk Nb sample, to their measured values. In the thin Nb 
layer, we set the effective correlation length and the average 
Fermi velocity to standard values, respectively £ s = 10 nm 
and v s = 4.4 x 10 5 m/s [24]. The dimensionless Ho exchange 
field ho, which determines the ferromagnetic coherence 
length £h , is an adjusted parameter. The calculated T c is 
shown in fig. 3 as a function of duo- To account for the 
measured magnetic dead layer, an horizontal shift ~ lnm 
has to be applied when comparing to the experimental 
data [17, 18]. We obtain, as in the experiments, multiple 
oscillations, with a similar T c decay and oscillations am- 
plitude, broadly reproducing the experimental shape of T c 
as a function of duo- Although the absolute value of T c 
differs (by ~ 1.5 K, see caption of fig. 3), such shifts have 
been previously seen in [21] due to the actual samples not 
quite conforming to the ideal assumptions of the theory. We 
find that the results are very sensitive to the value of d s , 
making it therefore remarkable that reasonable agreement is 
found without adjusting this parameter from its experimental 
value. The dn range of the T2 series is rather extensive 
and includes the small duo range where the T c oscillations 
are more pronounced. Of particular interest is the sharper 
minimum found near 10 nm. According to the theory this 
is due to the near emergence of a reentrance region in duo'- 
the system is very close to being non- superconducting at 
this thickness value, as found experimentally (see arrow in 
fig. 3). The adjusted values of the parameter ho used here 
give £ho — 3.2 nm for the Tl series and = 2.1 nm for 
the T2 series. These are both consistent with = 4.3 nm 
found experimentally [18]. Due to the ratio ^Ho/^ ~ 1, the 
oscillations in T c are governed by a combination of both the 
helical period A and the ferromagnetic coherence length £h q , 
leading to an intricate global behavior. While to provide a 
definitive test of the theory requires a knowledge of the exact 
magnetic structure, our primary aim is to demonstrate that 
our unusual results are at least consistent with a theory based 
on bulk Ho properties. 

R(T) dependence - We now consider the two separate tran- 
sitions appearing in the R(T) curves. To try to understand 
their origin, we have (i) applied a magnetic field up to 5 T 
and (ii) reduced the lateral dimensions of the films down to 
160 nm. As a function of the Ho thickness, we observe dif- 
ferent behaviors as in fig. l(a)-(c). If the transition is a single 
step at B = 0, the magnetic field can either act only as a 
pair-breaking mechanism, resulting in a progressive decrease 
of T c [fig. 1(a)], or it can split the transition in two resistance 
drops, whose distance in temperature increases with the field 
amplitude [fig. 1(b)]. If two transitions are already present at 
B = 0, they will grow apart with the magnetic field [fig. 1(c)]. 
In all cases, the amplitude of the two resistance drops is not 
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FIG. 4. (Color online) (a) T c ,io( H ) and T C M H ) for a Nb thin film 
(black lines, higher critical temperatures), and for two representative 
samples of the Tl series: duo — 5nm (red dots) and d,Ho = 6.2 nm 
(blue squares). Relative decrease with field of T c ,98 (b) and T c ,io (a) 
as a function of Ho thickness for B=0, 1,2,3,4,5 T from top to bottom. 
The scale is exactly the same in (b) and (c). 



affected by the magnetic field. When applying an out-of-plane 
magnetic field, the result is qualitatively the same, but the T c 
suppression is more than twice as strong. We have compared 
the in-plane magnetic field dependence of the two transitions 
with the field dependence of a control Nb layer. In fig. 4 (a) 
we trace T C5 io and T c $% for the control Nb layer and two rep- 
resentative bilayers from the Tl series. While T Cj 98 follows 
the dependence of the control Nb dependence with field, even 
though with a lower T c , T c ^ is suppressed on average 1.5 
times faster. Fig. 4 (b) and fig. 4 (c) show the relative decrease 
of T c with the magnetic field for the two transitions (T c 98 and 
T Cj io respectively) as a function of the Ho thickness. As al- 
ready observed in fig. 1, the magnetic field has different ef- 
fects over the two transitions: on average T c 98 is less affected 
by the field and its decrease is quite homogeneous over all the 
samples, while T C5 io decrease is larger and, most importantly, 
strongly dependent on the Ho thickness. 

We can thus identify two regions in our bilayers: in one, 
proximity superconductivity can be easily induced in the 
Ho layer, and the magnetic field affects similarly all the 
samples through the reduced T c of Nb; the other becomes 
superconducting at a substantially lower temperature, and the 
effect of the magnetic field is not limited to the Nb layer but 
also affects the magnetic structure for certain Ho thicknesses. 
The difference between the two transitions can be as high 
as 1.3 K, and is increased by the magnetic field, attaining 
in some samples 4 K. Since such two-phase behavior is not 
present in the Nb layer, we attribute its origin to the Ho layer 
and its magnetic structure. It can be explained if we assume 
the presence of magnetically-inactive regions which give rise 
to a higher temperature transition. Since the roughness is 
much smaller than the layer thicknesses [17], the Ho films 
are homogeneous, while the magnetic structure is not - either 
because of the presence of misoriented grains or because 
of domain structures. The size of the domains has been 
measured in epitaxially grown Ho to be ~ 50 fim [25], but 



we expect it to be smaller in our poly crystalline films. To 
explore the spatial dependence of the phase separation, we 
have drawn tracks on the film and made constrictions of fixed 
length L = 10 fim and variable width w in the 160 nm - 2 /xm 
range, by Focused Ion Beam etching. In the constrictions, 
the R(T) curves are qualitatively similar to those measured 
in the films, presenting two separate transitions at exactly the 
same critical temperatures as in the films, but with a varying, 
lower, step resistance (see [17]). This suggests that we can 
modify the proportion of the two phases without affecting 
their electronic properties when the device size approaches 
the characteristic size of the regions in which the phases are 
confined. 

Conclusions - In qualitative agreement with the theoretical 
modeling of inhomogeneous F-S junctions, our experiments 
show multiple complex oscillations of T c vs dn resulting 
from both the contributions of singlet and triplet Cooper 
pairs, and controlled by the helical magnetism and the ferro- 
magnetic coherence length in Ho. Even though the model is 
based on a monocrystalline structure, the agreement is sur- 
prisingly good, giving evidence of interesting and potentially 
important behaviour which will hopefully stimulate further 
experimental and theoretical research. The presence of a 
behavior which is already quasi reentrant in <i# , and which 
is in agreement with the theory, is particularly intriguing and 
will be addressed in future work. We furthermore observe 
that the R(T) curves show for certain Ho thicknesses two 
transitions, affected differently by the magnetic field. Further 
work is thus in progress to define the nature of the two phases 
in which the proximity effect acts so differently. 
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In this Supplementary Materials we provide the magnetic and structural characterization of our 
sputtered polycrystalline Nb/Ho/Cu trilayers. Moreover, we show supplementary T c vs. Ho thick- 
ness data in which we explore the high Ho thickness range. Finally, we show the effect of a lateral 
constriction on the films' R{T) curves. 
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I. HO MAGNETIC PROPERTIES 

The magnetization of the Nb-Ho films was determined from vibrating sample magnetometry at varying temperature 
for two Ho layer thicknesses, dn — 20 nm and dn — 100 nm [Fig. 1 (a)]. In all measurements the applied field was 
in-plane and along the long axis of the sample. The films become saturated at B w 4T. Extrapolation from the plot 
of the saturation magnetization as a function of Ho layer thickness reveals that the Ho has a magnetically dead-layer 
of ~ 1 ± 0.5 nm per Ho surface. The determined magnetisation is 2.7 x 10 6 A/m, lower than the theoretical value 
of 3.08 x 10 6 A/m, suggesting that magnetic inhomogeneity is present in our sputtered polycrystalline films. These 
measurements are in good agreement with the characterization of ref. [1]. 
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FIG. 1: (Color online) (a) Magnetic moment vs magnetic field at T=10 K for the two samples Nb(4)/Ho(20)/Au(10) and 
Nb(4)/Ho(100)/Au(10), where the Au layer protects the Ho from oxidation (the layer thicknesses are noted in nm units in 
the parentheses). The small jump at zero field is a measurement artifact, (b) Magnetic moment vs magnetic field for the 
Nb(4)/Ho(20)/Au(10) sample at T=7, 10,25,50 K. The slightly higher value of saturation moment in (a) is due to a better 
subtraction of the diamagnetic background. 
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II. SAMPLES CHARACTERIZATION 

Deposition rates and film thicknesses have been checked on pre-deposited step-edges and from low-angle X-ray 
reflectivity on Tl and T2 samples (Fig. 2); the latter gives an estimation of the film roughness of about 1.5 nm. 

Analysis of the integrated intensities of high-angle X-ray diffraction data taken from a thick bilayer film reveals 
that the Ho exhibits texturing with (002) as the primary phase, the c-axis being normal to the sample surface; a 
secondary (101) phase and a much weaker (100) tertiary phase are also present [1]. 

The Residual Resistance Ratio (RRR), the ratio of the resistance at 292K and 10K, was also measured to estimate 
the films quality. We obtained RRR~1.3, with a 20% variation from sample to sample. The samples of the two series 
Tl and T2 have similar RRR values, as expected from using the same fabrication method. 
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FIG. 2: (Color online) Low angle X-ray reflectivity for duo — 5nm of the Tl series (blue, higher curve). The theoretical 
curve (red, lower curve) is calculated for the multilayer Nb(14.4)/Ho(5)/Cu(17)/CuO(2.3) on a Si/Si02 substrate (the layer 
thicknesses are noted in nm units in the parentheses). The sample roughness of 1.7 nm is mainly due to the CuO layer, while 
the roughness at the Nb level is 1 nm. The disagreement at the lower angles may be due to little Nb oxidation at the Nb/SiCb 
interface. This, and the formation of a 1 nm dead layer, account for the lower Nb thickness found in the X-rays measurements. 
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III. RESISTANCE VS TEMPERATURE FOR HO SAMPLES 



Control Ho films, deposited over thin (4nm) Nb layers, show a monotonic R(T) curve in the experimental temper- 
ature range. Since neither the Ho films nor the Nb films present any heterogeneous behaviour, the double transition 
observed in the Nb/Ho bilayers can be attributed to the effect of the Ho magnetism on the Nb proximity effect. 
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FIG. 3: (Color online) Resistance vs. temperature for the Nb(4)/Ho(20)/Au(10) sample. The resistance drop at T=3 K is due 
to the superconducting transition of the thin Nb layer. Note that in the temperature measurement range 4-8 K, the resistance 
of the Ho layer is extremely stable. 



IV. LARGER HO-THICKNESS RANGE DATA 

To explore a larger range of Ho thicknesses, we have deposited Nb(15)/Ho(x)/Cu(20) trilayers with Ho thickness in 
the 4-32 nm range. The oscillations are here superimposed on a slow oscillation with a minimum around dn = 21 nm. 
This confirms the hypothesis that the slow decay seen in series Tl and T2 is but the initial part of a slow oscillation 
similar to that seen in homogeneous S-F bilayers. 
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FIG. 4: (Color online) Superconducting critical temperature T c ,5o vs. Holmium thickness d,Ho of a larger range series, fabricated 
in the same way as Tl and T2. The curves which delimit the shaded area are T C) io and T c ^8- 
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V. RESISTANCE VS TEMPERATURE FOR CONSTRICTIONS OF DIFFERENT WIDTH 

We show here the R(T) renormalized curves for a film in which constrictions have been created (see inset). While the 
amplitude of the two transitions is strongly and non- monotonously affected by the decrease of the lateral dimensions, 
we observe that the transitions temperatures correspond to those of the initial film. 




T(K) 

FIG. 5: (Color online) Renormalized resistance vs. temperature for duo — 7.75 nm of Tl. From top to bottom: film and 
constrictions with L — 2 jim and w — 1200 nm, w — 160 nm and w — 1000 nm. 
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